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Abstract 

We pose the question of what is the best generalization of the factorial and the 
binomial coefficient. We give several examples, derive their combinatorial properties, 
and demonstrate their interrelationships. 


Generalisation des Coefficients dn Binome 

On cherche ici a determiner est la meilleure generalisation possible des factorielles 
et des coefficients du binome. On s’interesse a plusieurs exemples, a leurs propriety 
combinatoires, et aux differentes relations qu’ils mettent en jeu. 

Dedicated to 
David and Maureen 
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Despite being so fundamental to combinatorics, several authors have noticed that ne is 
virtually unlimited in the choice of definition for the factorial —at least as far as umbral 
calculus is concerned. Indeed, one is presented with a bewildering number of alternatives 
each with its own notation. 

We present a new definition of the factorial which generalizes the usual one, and study the 
binomial coefficients it induces. They are blessed with a variety of combinatorial properties. 
However, what we are most interested is studying the interrelationship between this factorial 
and other famous ones. 
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1 INTRODUCTION 


1.1 The Roman Factorial 


We begin by presenting a generalization of the factorial n\ which makes sense for negative 
integral values of n as well as nonnegative called the Roman factorial [n]! after its inventor 
Steve Roman. As usual for n a nonnegative integer the factorial is given by the product 


[n"l ! = n! = 1 X 2 X 3 X ■ ■ ■ X n. 


However, for n a negative integer 


[n]! 


(_l)n+l 

(—n — 1)! 


Table 1: Roman Factorials [n]! 


n 

-6 

-5 

-4 

-3 

-2 

-1 

0 

1 

2 

3 

4 

5 

6 


n 


1 

120 

1 

24 

1 

6 

1 

2 

-1 

1 

1 

1 

2 

6 

24 

120 

720 


Proposition 1.1 (Knuth) For any integer n, 

[n]! [—n]! = (—l)”|n|. I 


More generally, for every real number a, let 

{ r(a + 1) when a is not a negative integer, and 

(—1)“ ^/{—a — 1)! when a is a negative integer 
where r(a) is the analytic Gamma function. 

Thus, for all a 

La]!/[«-!]!= La] (1) 

where Roman a is defined to be 

{ a for a 7 ^ 0 
1 for a = 0. 

Note that equation o and the condition [0] = 1 completely characterizes the Roman 
factorial of integers. 








1.2 The Roman Coefhcients 
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1.2 The Roman Coefficients 

These extensions of the notion of factorial leads to a corresponding generalization of the 
dehnition of Binomial Coefhcients. 


Definition 1.2 (Roman Coefhcients) For all real numbers a and b, define the Roman 
coefficient (read: “Roman a ehoose b”) to be 


a1 [a]! 

_b [b] \ [a — b]\ 


Table 2: Roman Coefficients, 


n 

k 


n\k 

-4 

-3 

-2 

-1 


2 

3 

4 

5 

6 

t ) 

- 1/840 

1/252 

- 1/56 

1/7 

1 6 

15 

20 

15 

6 


5 

- 1/504 

1/168 

- 1/42 

1/6 

1 5 

10 

10 

5 


1/6 

4 

- 1/280 

1/105 

- 1/30 

1/5 

1 4 

6 

4 

1 ( 

1/5 

- 1/30 

3 

- 1/140 

1/60 

- 1/20 

1/4 

1 3 

3 

1 r 

1/4 

- 1/20 

1/60 

2 

- 1/80 

1/30 

- 1/12 

1/3 

1 2 

1 r 

1/3 

- 1/12 

1/30 

- 1/60 

1 

- 1/20 

1/12 

- 1/6 

1/2 


1/2 

- 1/6 

1/12 

- 1/20 

1/30 

0 

- 1/4 

1/3 

- 1/2 

1 

11 1 

- 1/2 

1/3 

- 1/4 

1/5 

- 1/6 


-1 

1 

-1 

1 

1 -1 

1 

-1 

1 

-1 

1 


3 

-2 


-1 

1 -2 

3 

-4 

5 

-6 

7 


-3 


- 1/2 

- 1/2 

1 -3 

6 

-10 

15 

-21 

28 

-4 

1 f 

- 1/3 

- 1/6 

- 1/3 

1 -4 

10 

-20 

35 

-56 

84 

-5 

- 1/4 

- 1/12 

- 1/12 

- 1/4 

1 -5 

15 

-35 

70 

-126 

210 


When the two argnment are both integers, the relationship between the Roman coeffi¬ 
cients and the binomial coefficients is given by the following: 


Proposition 1.3 (The Six Regions) Let n and k be integers. Depending on what region 
of the Cartesian plane the point (n, k) is in, the following formulas apply: 


Region 1 lfn>k>0, then 


n 

.k 













Table 3: Region 1 


Table 4: Region 2 
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0 

1 

2 

3 

4 

5 

6 

7 

1 

7 

21 

35 

35 

21 

7 

1 

1 

6 

15 

20 

15 

6 

1 


1 

5 

10 

10 

5 

1 



1 

4 

6 

4 

1 




1 

3 

3 

1 





1 

2 

1 






1 

1 







1 









nV/c I 0 
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Region 2 If k > 0 > n, then 
Region 3 If 0 > n > k, then 




Table 5: Region 3 


n\k 

-6 -5 -4 -3 -2 -1 

-i 

-1 1-1 1-1 1 

-2 

5-4 3-2 1 

-3 

-10 6-3 1 

-4 

10 -4 1 

-5 

-5 1 

-6 

1 


Region 4 If k > n>0, then 


n 


-1 


= (-1) 
= (-1) 


n-\-k 


n+fc+1 


n — k\nj 
l(k-l' 


k\ n 


= (-l) 


n+k+l 


n + 1 


-1 


/ . ^ 

\n+lj 


-1 


i>o 


n-\-k 


A" 


X — k. 


x=0 


where the S{j,n) are the Stirling numbers of the second kind, and A is the forward 
difference operator Ap{x) = p(x + 1) — p{x). 

Region 5 lfn>0>k, then 


n 


= (-i)V 


I fn — k'' 


-1 


k\ n 


= (-l) 


— k — l' 


-1 


- 1 ) 


k+l 


— k — l'' 


k — n 
-1 


n 


n + 1 \ n+1 


n 


= —B{k — n,—k) 


A" 


X - n + k\x=o 




O to 00 
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Table 6: Region 4 



Table 7: Region 5 


n\k —4 

-3 

-2 

-1 

b -1/840 

1/252 

-1/56 

1/7 

5 -1/504 

1/168 

-1/42 

1/6 

4 -1/280 

1/105 

-1/30 

1/5 

3 -1/140 

1/60 

-1/20 

1/4 

2 -1/80 

1/30 

-1/12 

1/3 

1 -1/20 

1/12 

-1/6 

1/2 

0 -1/4 

1/3 

-1/2 

1 
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where the pair {n,n — k) lies in region 4 (defined above), and B{n,k) is the analytic 
Beta function. 


Region 6 Region 6: If 0 > k > n, then 


n 

.k 



1 f-n-l 


-1 


n + l\—k — Ij 

k — n — 1 


A 


k—n—l 


\k+l 


-n 


= (-l) 


fc+i 


x + n + l\x=o 
k — n — l 
k 


where the pair {k—n—l, —n—1) lies in region 4 (defined above), and the pair {k—n—l, k) 
lies in region 5 (defined above). I 


Table 8: Region 6 


n\k 

-6 

-5 

-4 

-3 

-2 

-1 

-2 






-1 

-3 





- 1/2 

- 1/2 

-4 




- 1/3 

- 1/6 

- 1/3 

-5 



- 1/4 

- 1/12 

- 1/12 

- 1/4 

—6 


- 1/5 

- 1/20 

- 1/30 

- 1/20 

- 1/5 

-7 

- 1/6 

- 1/30 

- 1/60 

- 1/60 

- 1/30 

- 1/6 


Note that in regions 1, 2, and 3, the Roman coefficients equal binomial coefficients up 
to a permutation and a change of sign. In regions 4, 5, and 6, the Roman coefficients 
are expressed simply in terms of the reciprocals of the binomial coefficients. Furthermore, 
regions 4, 5, and 6 are identical up to permutation and change of sign. Thus, all of the 
Roman coefficients are related in a simple way to those in the hrst quadrant (regions 1 and 
4). In particular, the Roman coefficients always equal integers or the reciprocals of integers. 

1.3 Properties of Roman Coefficients 


Several binomial coefficient identities extend to Roman coefficients. 
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Proposition 1.4 (Complementation Rule) For all real numbers a and b, 


a' 


a 

_b 


.a — b 


Proposition 1.5 (Iterative Rule) For all real numbers a, b, and c; 


a' 


h 


a' 


a 

— c 

b 


c 


.c 


b 

— c 


Proposition 1.6 (Pascal’s Recursion) If a and b are distinct and nonzero real numbers, 
then we have 

a1 a — 11 a — 1 
.k k ^[k-1 


Proof: Since under these conditions [a] = a, [b~\ = b, and [a — b] = a — b, 


a — 1 
b 


a — 1 
b-1 


[a - 1]! 


+ 


[a - 1]! 


[a — b — 1] \ [b]\ [a — 6] ! [6 — 1 ] ! 


= L«1 


[a — b] \ [b ]! ^ 
[a - 1]! \ 


[a — b] \ [b ]! _ 


\ [a — 6]! [b ]! 
La]! 


[a — b] \ [b ]! 
a' 


Corollary 1.7 Ifr is a nonnegative integer and the pairs of integers (n, k), {n + r, k), 
(n, k + 1), and {n + r + 1, k + 1) all lie in the same region (as defined in Theorem n) , we 
have 


n+r 

E 


m=n 


m' 


n + r + 1" 


n 

k 


- k + 1 


k + 1 
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Proof: Induction on r. I 


Contrast this corollary with this classical result involving binomial coefficients in which 
for n > fc > 0, 


E 




Analogous results hold more generally for real numbers. 

If we adopt Iverson’s notation for the moment writing logical expressions in parenthesis 
to mean 1 if true and 0 if false, then in the discrete case we have the following proposition: 


Proposition 1.8 (Knuth’s Rotation/Reflection Law) For any integers n and k, 


- l ) fc +( A :> 0 ) 


—n 

^k-l 


(- 1 ) 


n+(n>0) 


-k 
n — 1 


Proof: By Proposition |1.1| , we have 


= (-l) 


n+fc+(n<0)+(fc<0) 


-k-1 
—n — 1 


Proposition 1.9 (Roman’s Identity) For all integers n and k, 


(- 1 ) 


n-\-k 


\n — k\ 


Proof: Proposition 0- ■ 
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1.4 Generalizations of the Roman Coefficients 

The Roman coefficients defined earlier were very useful. However, there are several other gen¬ 
eralizations of binomial coefficients. For example, recall the classical definition of extended 
binomial coefficients. 


Definition 1.10 (Classical Extended Binomial Coefficient) Given a field element x G K 
in a field K of characteristic zero, and a nonnegative integer k, define the binomial coefficient 
“x choose k” to be: 

(l) = 

where {x)k denotes the lower factorial of x of degree k 


{,x)k 


k-l 


n(^-o 

i=0 


n(^-o ^ 

. i=k 


x{x — T) ■ ■ ■ {x — k + T) for k>D, and 
1/{x + l){x + 2) ■ ■ ■ {x — k) for k < 0. 


What is the relationship between the Roman coefficients and the other generalizations 
of binomial coefficients? To fully answer this question, we must generalize our notion of 
harmonic factorial. 


Knuth Coefficients 


Adopt the following convention independently discovered by Donald Knuth. 


Definition 1.11 (Knuth Factorial) Define [a]^ for a a real number to be the most 
significant term ofT{a + l + e) where e is an infinitesimal in from the field of surreal numbers 
(a non-Euclidean field which contains the real numbers). 
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Thus, for a real, 


L«i: 


r(a + 1) when a is not a negative integer, and 

(—l)““^ci;/(—a — 1)! when a is a negative integer. 


( 2 ) 


where a; = 1/e. This choice of factorial whoold have led to “tags” of e or a; in appropriate 
places in results of this paper. 


For instance, again for a real. 




a if a 7 ^ 0, and 
e if a = 0. 


This is perhaps more natural since then [a] ^ only differs from a by at most an infinitesimal. 

If we adopt equation (^) as our definition where e can be any arbitrary constant, then 
the Roman factorial can be seen as a special case of the Knuth factorial where e = 1. That 
is, [a]! = [al^. Thus, the motivation for our notation. 

Let us proceed to generalize the Roman coefficients. 


Definition 1.12 (Knuth Coefficient) For all a and b, define the Knuth coefficient 
by the fraction 


n: 


[K [a-byfi 


Clearly, “ 


Let us calculate 



for each of the six regions mentioned in Theorem |1.3| . 


Proposition 1.13 (The Six Regions) Let e be a nonzero complex number or surreal 
number, and and n, k he integers. Depending on what region of the Cartesian plane the pair 
(n, k) is in, the following formulas apply: 

Region 1 lfn>k>0, then 

n 
.k 
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1 INTRODUCTION 


Region 2 If k > 0 > n, then 


= (-l) 



Region 3 IfO>n>k, then 


n 

.k 


(- 1 ) 


n+fc 


'-k-C 

, n — k , 


n 

.k 


Region 4 If k > n > 0, then 


n' 


n' 

.k 


.k 


Region 5 lfn>0>k, then 


n 

e (_l)n+fcg 

n 

.k 

•” TO 



Region 6 If 0 > k > n, then 



Gamma-CoefRcients 


A limiting case of the Knuth coefficient is of special interest. 


Definition 1.14 (Gamma-Coefficient) Let n and k he arbitrary integers. Define the 
Gamma- Coeffieient 


n 

0 

n 

= lim 

.k 

^ e —>^0 

.k 


= lim 


r(n + 1 + e) 


0 r(A; -h 1 e)r(n -k + l + e) 
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Table 9: Gamma-Coefficient 


n 

k 


0 


n\k 

-4 

-3 

-2 

-1 

0 

1 

2 

3 

4 

5 

6 

b 

0 

0 

0 

0 

1 

6 

15 

20 

15 

6 

1 1 

5 

0 

0 

0 

0 

1 

5 

10 

10 

5 

1 

0 

4 

0 

0 

0 

0 

1 

4 

6 

4 

1 [ 

0 

0 

3 

0 

0 

0 

0 

1 

3 

3 

1 

0 

0 

0 

2 

0 

0 

0 

0 

1 

2 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

-i 

-1 

1 

-1 

1 

1 

-1 

1 

-1 

1 

-1 

1 

-2 

3 

-2 

1 

0 

1 

-2 

3 

-4 

5 

-6 

7 

-3 

-3 

1 

0 

0 

1 

-3 

6 

-10 

15 

-21 

28 

-4 

1 

0 

0 

0 

1 

-4 

10 

-20 

35 

-56 

84 

-5 

0 

0 

0 

0 

1 

-5 

15 

-35 

70 

-126 

210 


Note however that 


Gamma-Coefficient ^ for a and b real. 


l0 


-1 

[ 1/2 


diverges as e tends to zero, so it is impossible to dehne a 


In regions 1, 2, and 3, the Gamma-Coefficients are equal to the Roman coefficients. 
In regions 4, 5, and 6, the Gamma-Coefficients are identically zero whereas the Roman 
coefficients are never zero. Nevertheless, one should note that even when the Classical 
binomial coefficient and the Roman coefficient differ, the difference is at most one. 


Also, notice that the Gamma-Coefficients are always integers. In particular, for k > 0 
{i.e.: regions 1, 2, and 3), the Gamma-Coefficients agree with the classical extended binomial 
coefficients. 


The identities mentioned in §H generalize to Gamma-Coefficients. However, we defer 
any discussion of the combinatorial signihcance to [0. 


a 

a—b 


Proposition 1.15 (Complementation Rule) For all real numbers a, b, and e, 
In particular, for all integers n and k, 


n 

0 

n 

_ k 

• 

n—k 
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1 INTRODUCTION 


Proposition 1.16 (Iterative Rule) For all real numbers a, b, c, and e, 


a' 

€ 

h 

e 

a' 

e 

a — c' 

_b 

• 

c 

• 

.c 

• 

_h — c 


In partieular, for all integers m, n, and k, 


m' 

0 

n' 

€ 

m' 

0 

m 

-k 

. n 

• 

k 

• 

.k 

• 

n 

-k 


Proposition 1.17 (Pascal’s Recursion) 

1. Let a and b be distinct nonzero real numbers, and let e be a nonzero complex number. 
Then 


aR 


a-lV 


a-lV 


b 


b 


b-1 


2. For all n and k, 





n — iR 
— 1 • 


unless n = k = 0. 


Nevertheless, 


= 1 whereas 


-1 

-1 


l0 


+ 


-1 


l0 


1 + 1 = 2 . 


Other Factorials 


Actually as noted by Ueno |T^ and Roman |^, |^, any choice of [a]! could be used for 
computations involving an umbral calculus. The only restrictions are that [0]! must equal 
one, and for the so called continuous iterated logarithmic algebra of p, the function a ^ [a]! 
must be continuous. 

For example, if we chose [a]! = 1 as in []^, then we have the theory of convolution 
sequences. 


Whereas, if for n an integer, we set as in |12 


n = 


q 


+1 _ 1 


g - 1 ’ 

then we achieve a g-analog of the “[u]-Logarithmic theory.’ 
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Multinomial Coefficients 

Recall the usual definition of a multinomial coefficient. 

Definition 1.18 (Classical Multinomial Coefficient) Letn be a nonnegative integer, and 
let 13 be a vector with finite support of nonnegative integers. Then define the multinomial 
coefficient n choose [3 to be 

if \f^\ = n, and 
0 otherwise. 



Note that Qj is the number of ordered partitions of type (3 of a given n-set. 

By analogy, for all reals a, and all real vectors (3 with finite support, define the multinomial 
Roman coefficient a choose j3 to be 


a 


/3 


wifriL&i!) 

0 


if \j3\ = a, and 
otherwise. 


Define the multinomial Knuth coefficients and Gamma-Coefficients similarly. The multino¬ 
mial Gamma-coefficients are well defined since they would only diverge if some denominator 
had an excess of factors of e. However, that could only happen if n < 0 and fcj > 0 for all i, 
but in that case, n ^ Y3l=i ki, so the multinomial e-coefficient, y is zero by definition. 
Contradiction! Thus, the Gamma-coefficients are well defined. 


In terms of multinomial coefficients. Proposition |T]^ becomes 


n' 


k 


n 


n' 

n — r' 

.k 


r 


n — k,k — r,r 


_r 

_k — r 


and Proposition |1.16| becomes 


n' 

0 

k' 

0 

n 

0 

n' 

0 

n — r' 

k 

• 

r 

• 

n — k,k — r,r 

• 

_r 

• 

_k — r 


More generally, we have the following theorem. 
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Proposition 1.19 (Iterative Rule) Let be a finite sequence of integers with sum 

n. Then 


n 

{k^)U 

n 


J L, 

l^i=l 

kr. 


m=2 L ™- 

j 


n 

{ki)l=i 


n 0 


n 

m=2 

3 

n 

m=2 

j 

n 

m=2 


TZih 

k 

rvm 

YZih 

h 


and 


As opposed to ordinary Roman coefficients, these multinomial Roman coefficients are not 
always integers or reciprocals of integers—even when all of the arguments are integers. For 
example, 22^-1 = §• 


l0 


liWl 


IS 


However, the multinomial Gamma coefficients are always integers, for if 
nonzero, then we are in one of the following two cases. Either n > 0, and ki > 0 for 
all i, or n < 0 and there is a unique i such that fcj < 0. In the hrst case, these are ordinary 
multinomial coefficients. It suffices to consider the other case. Thus, n < 0. Without loss of 
generality, let fci < 0. Now, 


n 

iki)U 


n 0 


= (-l) 


n+fci 


-n 


—ki — 1 


n 0 


, kj 


= (-l) 


n+fci 


-ki-i 

-n-l,k2w,kj^ 


where —ki — 1, —n — 1, /c 2 , k^, - ■ ■, fcj-i, kj > 0. Hence, all the nonzero multinomial Gamma 
coefficients are (up to sign) ordinary multinomial coefficients, and thus integers. 


1.5 Resistance of the n-cube 

Via the Gamma-coefficients and the theory of sets with a negative number of elements 0, we 
have a simple combinatorial interpretations for the Roman coefficients in regions 1, 2, and 3. 
However, what is the signihcance of the Roman coefficients in regions 4, 5, and 6? In these 
regions, the Roman coefficients are the reciprocals of integers, so they do not enumerate any 
set. However, the following application illustrates their combinatorial signihcance. 
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Figure 1: The n-cube 


Proposition 1.20 Consider an n-cube in which each edge is represented by a wire of 
resistance IQ (one Ohm). The resistance between two opposing vertices of the cube is 

n 

i=l 


Table 10; Resistance of the n-cube 


n 

0 12 3 

4 

5 

6 

7 

Rn 

0115/6 

2/3 

8/15 

13/30 

151/340 


Proof: The cube is isomorphic to the Hasse diagram of the boolean lattice of subsets of 
{1,2,..., n}. Without loss of generality, the two opposing vertices are 0, and (1,2,..., n}. 
To compute the resistance, connect these two vertices to a IV battery. The resulting current 
(in Amperes) is equal to the resistance (in Ohms). 

By symmetry, each vertex on level i of the lattice has the same potential. Hence, we can 
consider each level as a single node without effecting the resistance. Any two adjacent levels 
i and i + 1 are connected by (n — i) ff) edges. Thus, the resistance between levels i and i + 1 
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is or in the notation of Roman coefficients, the resistance between levels i 

n—i^ ’ 


i + 1 is (— 


Q = - 


—n—1 

1-i 


n. 


The total resistance Rn is the sum of the resistances between the adjacent levels, 


-1 




—n — 1 
i 


n. 


By Theorem 


2-Rn — Rn-1 H- 

n 

We conclude by induction noting that Rq = 0 and Ri = 111. 


Note that as n tends towards inhnity, Rn tends towards zero as 2In. 


i and 
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